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ABSTRACT 

It is shown that the influences of the thermal and tidal effects on Mercury's 
libration are in an equilibrium condition with the periods of rotation and revolu- 
tion of Mercury locked in the 3:2 resonant state. 
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THERMAL AND TIDAL EFFECT ON THE 
LIBRATION OF MERCURY 


INTRODUCTION 

Liu (Liu, 1968a; Liu, 1968b) has shown that the trapping of Mercury's rota- 
tional period into a 3:2 resonance lock with its orbital period was originally 
affected by the thermal contraction of the figure of Mercury during solidifica- 
tion. In the present paper the analysis of the contribution of the two thermal 
bulges to the dynamic stabilization of the planet's libration is given. Attention 
is focused on the balance of the influence of the tidal and thermal effect on the 
libration of Mercury after solidification. 


BASIC EQUATIONS 

After solidification, the thermal contraction of the figure of Mercury from 
loss of heat must be exceedingly small and the thermal bulges considered by 
Liu (Liu, 1968a; Liu 1968b) can grow because the apparent circulational motion 
at successive perihelia has been converted to a librational motion. The varia- 
tion with time of the fractional difference [b ( t ) - A ( t in Mercury's equa- 
torial moments of inertia is 


B (t) ~ A (t) _ 3(2 ~ 5/3 3 + 3/3 5 ) 

C 100(l +/3+/? 2 ) 


(1) 


in which a is the coefficient of linear thermal expansion, AT the difference in 
surface temperature between the regions around the perihelion and the aphelion 
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axes of Mercury, and /3 ~ (R m -y) / /R m where R ra is Mercury's radius and y the 
depth of solar heating. The value of \ has been estimated as A. > 5 x 10“ 5 and 
the rate of increasing of ^B ( t ) - A ( t is about 10“ 2 1 sec. ' 1 (Liu and 
O'Keefe, 1965; Liu, 1968a; Liu, 1968b). 

In considering the bodily tidal torque we have recourse to the estimation by 
Jeffreys. (Jeffreys, 1959). 

1 18 M s 2 R m 6 (l +e cos f) 6 

N = §0 *? VQP ^ ‘ ' a 6(!- e 2)6 ‘C*Sin(20 (2) 

where G is the gravitational constant, p density, M s mass of the Sim, M m mass 
of Mercury, f the true anomaly, a the semimajor axis, e the orbital eccen- 
tricity and e the phase lag of the conventional equilibrium tides. 

The orientation of Mercury relative to the Sun, <£, is then governed by (Liu 
and O'Keefe, 1965) 


d 2 cp 2e Sin f /dp \ 3(A. +&V) _ N_ (1 ~e 2 ) 3 (3) 

df2 1 + eoosf ^dF + 1 J + 2(l+eoosf) -Sin 2P = " n 2 C ' ( 1 + e cos f ) 4 

in which 


A K 


3(2 ~5/3 3 +3/3 5 ) 

100(l +/3 +,3 2 ) 


and n is the mean orbital angular velocity. 
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ANALYSIS 

Let us choose the time r reckoned from perihelion and related to the 
orbital period divided by 2 n 


r ( f ) = 2 tan 


- 1 




/2 f e(l - e 2 ) l/2 Sin f 


tan "2 


1 +o cos f 


( 4 ) 


As a new unknown function, we take the angle of rotation ¥ between the axis 
of the two thermal bulges and the radius vector of the perihelion 


¥ = f + 0 


(5) 


Equation (3) takes the form 


d 2 <y 3(\ + A\) (ltecos f) 3 


dr 2 2 (l-e 2 )3 *' n 2 C 


Sin 2(» ~ f) = “ 

The solution of Equation (6) may be sought in the forrr. 


N 


( 6 ) 


= fir + 8 

where H is a constant and S is an unknown function. The resonances occur at 
2ft = m if m is an integer and 8 is, then, the angle of libration. To obtain an 
approximate solution we average it over a period of 2n. For m = 3, Equation 
(6) becomes 


d 2 S 
dr 2 


2 • 1 

+ — s— 


(K +A\) Sin(2S) 


10" 9 Sinf2e) 


(7) 
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In the derivation of Equation (7), the following values were adopted: 


e = 0.2 


G = 6.7 x 10" 8 dyn. Cm 2 • g~ 2 
p = 5.0 g • Cm -3 


n = 1.2 x 10“ 6 rad. • Sec' 


M./M = 6.0 x 10 6 


fi /a = 4.1 x 10" 


The average process also included the effect of the sign of the phase lag which 
depends on the sign of dtp/d r. In Equation (7) the value of K + AA. is about 
5 x 10" 5 and the tidal torque is much less than the maximum restoring thermal 
torque. Therefore, the rotation of Mercury at the 3:2 resonance state is seen 
to be stable. 

The first integral of Equation (7) is 


1 /d 8\ 2 2. 1 

J \SF) - - 4 - (K +AX) Cos( 2 S) 


where E n is a constant. 


Differentiating Equation ( 8 ) with respect to t , we obtain 


1 d (db\ 2 

2 3t \dt j 


2 • In - 5 . _ d Cos( 

—4— (^+A\) — dT 


2 • In 2 d(A\) 

= —4 W"-Cos(2 5 ) 


10“ 9 n 2 Sin( 2e ) ^ 
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Hence, 


dE dE (thermal effect) dE (tidal effect) ^ 

at = — 1 1 — ! ar l — at < 10 > 


in which 


dE (thermal effect) 

3t 


4 - n 2 • C • Cos( 2i > ) 


d(A\) 

•~3t~ 


and 


( 11 ) 


dE (tidal effect) „ „ 

= - 10” 9 n 2 • 

For small librational angle b and d(A\)/dt = 10" 21 
tion (11) is 


C • Sin( 2e ) -gt- (12) 

Sec" 1 , the result of Equa- 


dE (thermal effect) 
dt 


0 ( 10 10 erg. • Sec -1 ) 


For Sin(2e) - 5 x 10" 3 and dS/dt < 10” 3 n, the result of Equation (12) is 


dE (tidal effect) 
dt 


0 ( 10 10 erg. • Sec' 1 ) 


It is seen that the thermal expansion effect works against the tidal effect on 
Mercury's libration at a rate of 10 1 0 erg. Sec” 1 . This is of the same order of 
magnitude as the rate of the librational dissipation for the bodily tidal friction 
of Mercury. 
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CONCLUDING REMARKS 

We "nave shown that the influences of the i hermal and tidal effect on 
Mercury's libration are in an equilibrium condition with the periods of rotation 
and revolution of Mercury locked in the 3:2 state. With regard to the interaction 
of the thermal effect, tidal friction and gravitation, we conclude that the growth 
rav3 of thenr.nl bulges on Mercury's surface contributes dynamically to the 
stabilization ot the planet's libration. 
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